Phase Description of Stochastic Oscillations 
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We introduce an invariant phase description of stociiastic oscillations by generalizing the concept 
of standard isophases. The average isophases are constructed as sections in the state space, having 
a constant mean first return time. The approach allows to obtain a global phase variable of noisy 
oscillations, even in the cases where the phase is ill-defined in the deterministic limit. A simple 
numerical method for finding the isophases is illustrated for noise-induced switching between two 
coexisting limit cycles, and for noise-induced oscillation in an excitable system. We also discuss how 
to determine the isophases for experimentally observed irregular oscillations, providing a basis for 
a refined phase description of observed oscillatory dynamics. 
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Phase reduction is the basic tool in the characteri- 
zation of self-sustained, autonomous oscillators. With 
a reasonably defined phase variable, one obtains a one- 
dimensional representation of the oscillator, allowing to 
describe important aspects of its dynamics, such as reg- 
ularity, sensitivity to forcing and noise, etc [H-Q- Fur- 
thermore, the concept of phase reduction is substantial 
for the data analysis of oscillatory processes in physics, 
chemistry, biology, and technical applications, where var- 
ious approaches exist for extracting phases from oscilla- 
tory time series 0-13 • 

To understand many properties of oscillating systems, 
such as their phase resetting and synchronizability, it is 
important to define the phases not only for the purely 
periodic motion, but for the whole state space. In the 
theory of deterministic oscillations this is done via so- 
called isochrones 9], by attributing to a state with an 
arbitrary amplitude the phase of a point on the limit cy- 
cle, to which this state asymptotically converges. In this 
letter we generalize this concept to irregular, noisy oscil- 
lators. The main idea is based on the definition of the 
isophases by virtue of the mean first passage time con- 
cept. We will first apply our method to noise-perturbed 
deterministic oscillators for which the isophases can be 
compared with deterministic isochrones. Furthermore, 
we will consider examples for which the isochrones and 
even the oscillations themselves disappear in the deter- 
ministic limit. The method will be also applied to noise- 
perturbed chaotic oscillations. Finally, we will demon- 
strate its applicability to experimentally observed time 
series. 

We start by reminding the standard definition of 
isophases (which in this case are also isochrones) in de- 
terministic systems with a stable limit cycle Xo(t) — 
Xo{t+T) having period T. First, one defines the phase on 
the fimit cycle v(xo). Then, being observed stroboscopi- 
cally with time interval T, all the points x that converge 
to a particular point on the cycle Xq having the phase 
(/3(xo). These points form a Poincare surface of section 
J((/3(xo)) for the trajectories of the dynamical system, 
with the special property that the return time to this 



surface equals T for all points on it. Thus, to find an 
isophase surface is equivalent to find a Poincare surface 
of section with the constant return time T. 

For a noisy system we define the isophase surface J 
as a Poincare surface of section, for which the mean first 
return time J — > J, after performing one full oscillation, 
is a constant T, having a meaning of the average oscil- 
lation period. In order for isophases to be well-defined, 
oscillations have to be well-defined as well: amplitudes 
should always remain positive (so that one can reliably 
recognize a "full oscillation"). Otherwise the concept of 
phase is not well-defined for a random process, and the 
isophases are not meaningful. 

Analytical calculations of the mean first return time 
(MFRT) is a complex problem in dimensions larger 
than one, therefore below we apply a simple numer- 
ical algorithm for construction of the isophases: an 
initial Poincare section is iteratively altered until all 
mean return times are approximately equal. In two- 
dimensional systems for which isophases are lines, we 
represent Poincare sections by a linear interpolation in 
between a set of knots. For each knot xj, the average 
return time Tj is computed via the Monte Carlo simu- 
lation. According to the mismatch of Tj and the mean 
period (T), the knot Xj is advanced or retarded. The 
procedure is repeated with all knots, until it converges 
and all return times Tj are nearly equal to (T) . 

Before proceeding to different applications of this pro- 
cedure, we discuss the importance of knowing isophases 
for noisy oscillations. The first important application 
is that of phase resetting. Phase resetting curve deter- 
mines how an oscillator responds to an external kick, this 
response determines synchronization properties of oscil- 
lator 13, HI- For deterministic oscillators the phase re- 
sponse curve is determined just from the isochrone to 
which the kick shifts the state of the system from the 
limit cycle. For irregular oscillators the proper definition 
of the phase response curve is based on the first passage 
time 121, so to determine it one has to find to which 



isophase, as defined above, the system is shifted by the 
kick (note that in the limit of small noise, perturbative 
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approaches for the phase dynamics do the job 13l. Il4|). 
The second apphcation is in the analysis of experimental 
data of coupled oscillators (cf. There one needs to 

determine the phase dynamics from the time series, this 
task is relatively simple to accomplish if the variations of 
the amplitudes are very small so that the definition of a 
phase-like variable along the observed limit cycle is un- 
ambiguous. However, in the presence of large irregular 
amplitude variations, the phase characterization of the 
oscillations is not unique (cf. Fig. |6] below) . Proceeding 
according to the given above definition of the isophases 
as the lines on the two-dimensional embedding plane, for 
which the mean return times do not depend on the am- 
plitudes, allows us to get rid of the ambiguity and to 
determine the phase in a consistent way. 

We stress here that in our definition of the isophases 
we do not assume Markovian property of the process, if 
the dynamics is non-Markovian, then the definition of 
the MFRT includes averaging over the "prehistory" or 
hidden variables as well. To illustrate this we consider 
as the first example a simple Stuart-Landau oscillator 
(variables r, 9) perturbed by an Ornstein-Uhlenbeck noise 
(OUN) C(i): 



r = r(l - r^) + (TrC{t) 

ic = -c+ v^m , 



K(r 



I) 



(1) 



where £^{t) denotes a (5-correlated white noise, 7 is the 
correlation time of the OUN, cu is the frequency of the 
noise-free limit cycle, and k is a nonisochronicity param- 
eter. In the state space r, 6, C the process is Markovian, 
but on the two-dimensional plane r, 9 it is not. Nev- 
ertheless, by the method described we obtain numer- 
ical isophase for which the MFRT is nearly constant 
(Fig. [Ij. This isophase can be obtained also from the 
following analytic approximation. First, we introduce a 
"corrected" phase variable ip = 6 — nlnr which obeys 
= Lo + crK(7C — Averaging this expression and 

identifying lo — ip where ip is the correct uniformly rotat- 
ing phase, we obtain Lp = tp — aKj(^. In this expression we 
have to account for correlations of C and r, to obtain the 
isophases on the plane {0,r). Assuming that r follows 
C(t) adiabatically, we obtain crC « — 1, what leads to 
the following expression for the isophases 



(p = 9 — tilnr — Kj (r^ — l) 



(2) 



Isophase following from this formula is compared with 
numerical one in Fig. [T] Interestingly, the noise-induced 
correction (last term in ([2])) does not contain parameter 
a, but the range where this correction is valid |f ^ 1| ^ c 
shrinks with the noise amplitude. 

While in the simplest example above the effect of noise 
is in the correction of the deterministic isochrones only, 
we consider now a situation where local isophases of dif- 
ferent periodic motions are "mixed" by noise resulting in 
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FIG. 1. (color online) The average isophase for noise-driven 
oscillator [Eq. [T] with oi = 1, k = 1, cr = 0.15, 7 = 1 in carte- 
sian coordinates. The inset shows the reduction of the vari- 
ations of the MFRT compared to an arbitrary cross-section. 
The numerically-derived isophase shows minor differences to 
the analytic approximation [Eq. [5] . 



new, global isophases. To this end we analyze the follow- 
ing model of two coexisting stable limit cycles, driven by 
white noise: 



r = r{l- r)(3 - r)(c ~ r) + a£,{t) 
9 = uj + S{r-2)-{l-r){3-r) . 



(3) 



Without noise, the system shows two limit cycles r/ = 
l,rjj — 3 (which have the same frequency if S — 0), 
separated by an unstable cycle at r = c. Each of the 
stable cycles has its own isophases, which meet singu- 
larly (as infinitely rotating spirals) at the basin bound- 
ary r — c. With noise, trajectory switches between the 
basins, so that combined mixed-mode oscillations involv- 
ing both cycles occur. By applying our method, we find 
the isophases of these oscillations in the whole range of 
the amplitudes, as shown in Fig. [2] While for small noise 
amplitude a residue of the singularity at the basin bound- 
ary is clearly seen, for a strong noise the isophases are 
rather smooth curves. 

Another example where otherwise singular isophases 
are smeared by noise is that of chaotic oscillations. Many 
chaotic attractors allow a representation in terms of am- 
plitudes and phases d, [H, [l^ , but because the phase 
generally performs a chaos-induced diffusion, isophases 
in the strict sense do not exist. Recently, description of 
chaotic oscillations in terms of approximate isophases has 
been suggested 17|. With noise, the return times to a 



Poincare surface of a strange attractor can be defined in 
the averaged sense only, and in this respect there is no 
difference between chaotic and regular deterministic os- 
cillators. Thus, the procedure of finding isophases based 
on the constancy of the MFRTs can be applied to chaotic 
systems as well, as is illustrated in Fig. |3] for the Roessler 
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FIG. 2. (color online) Isophases of bistable oscillations mixed 
by noise, in model ((3]) with oj = 3, c = 1.8, 5 = and different 
noise intensities. Solid red curves are stable cycles, the dashed 
curve is the unstable cycle representing the basin boundary. 

model 

i = -y - z + (T£,i{t) , 

y = x + 0.16y + a6W , (4) 
i = 0.2 + z{x - 10) , 

with uncorrelated Gaussian white noises in x and y com- 
ponents. 
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FIG. 3. (color online) Isophases of the noise-driven Roessler 
chaotic system Q, for different noise intensities u. Smaller 
noise leads to less smooth curves. Grey dots show the deter- 
ministic Roessler attractor. 

Our final example are noise-induced oscillations in an 
excitable system, which without noise has just a stable 
steady state, so deterministic isophases do not exist in 
any sense. With noise, such a system demonstrates oscil- 
lations which may be quite regular in the case of coher- 
ence resonance [18|. To build the model, we modify the 
noisy Stuart-Landau oscillator, with y- polarized noise, 
to perform noise-induced oscillations: 

r = r(l - r^) + err cos 9 £_{t) , 

9 (5) 
6 ^ uj + rcos9 ~ Kr + as'mO ^{t) . 



For < — K < 1 there is a stable state at tq = 1, = 
TT — arccos(aj — k) and an unstable state at ri = l,6i — 
TT + arccos(a; — k) (at — k = 1 they give rise to a periodic 
oscillation via a SNIPER bifurcation). Noise (cr 7^ 0) 
excites the state ro, 6*0 beyond ri, 0i and produces noise- 
induced oscillations (Fig. H]) . For strong excitability and 
small noise, the phase is well- defined, and the isophases 
can be introduced as curves with constant MFRTs. We 
show ten isophases in Fig. [5] (examples Fig. [Hand Fig. [2] 
have been rotationally symmetric, so one drawing of one 
isophase was sufficient, here the rotational symmetry is 
broken) . The effect of the noise intensity on the isophases 
is maximal at ~ tt, i.e. in the region of excitability 
where the oscillations spend most of the time; in the 
"deterministic" region 16*1 < tt/2 the isochrones are less 
sensitive to noise. 
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FIG. 4. At cj = 1.99, K = I, and a = 0.6, system Q shows 
seemingly self-sustained oscillations (with some rather noisy 
patches as well), where indeed oscillations are noise- induced. 
Top panel: variable x{t) = rcosO, bottom panel: 0{t). 
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FIG. 5. (color online). Average isophases of Eq. ((5]) at 
oj = 1.99 and k = 1 vary with noise intensity a as indicated. 
Larger noise intensity makes the system less isochronous, let- 
ting average isophases show stronger curvature. (Background 
trajectory of noise-induced oscillations (grey line) corresponds 
to the strong noise case a = 0.6.) 

A practical definition of the isophases for which the 
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MFRT is constant, is straightforward for numerical mod- 
els of irregular oscillators as illustrated above, but it can 
be used for experimentally observed signals as well. For 
this purpose one needs a two-dimensional embedding of 
observed oscillations, which can be, e.g., achieved by us- 
ing the Hilbert transform of the signal as the second vari- 
able. In Fig. |6] we present such a representation of mea- 
surements of human respiration, taken from the Phys- 
ionet database [l^l- One can see that the oscillations 
have a large amplitude variability, and defining the phase 
has a large degree of ambiguity - contrary to the situa- 
tions with a nearly constant amplitude, where a similar 
embedding on the signal vs its Hilbert transform plane 
results in a very narrow band of trajectories. The initial 
phase-like variables and the isophases resulting from the 
iterative procedure as described above are presented in 
Fig. [6] Application of the calculated isophases to deter- 
mining the phase dynamics of the observed signals gives 
the mostly uniformly rotating phase, maximally uncor- 
related from the amplitude variations. 
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FIG. 6. (color online). Left panel: Arbitrary phaselike vari- 
ables (black) and isophases (red) of the experimental data of 
a respiration signal taken from the public database "Phys- 
ionet". Right panel; for isophases the return times are inde- 
pendent of the radial variable (cross-correlation vanishes). 

Summarizing, we have introduced for irregular oscilla- 
tions a concept of average isophases, based on the con- 
stancy of the mean first return times. By applying a sim- 
ple procedure, we determined these isophases in a unified 
way for different classes of noisy oscillators: (i) noise- 
perturbed periodic oscillators, which posses isophases 
also in the noise-free case; (ii) multistable oscillators 
which in the noise-free case posses different singular 
isophases, but the latter become well-defined when dif- 
ferent modes merge due to noise; (iii) chaotic attractors 
where in the purely deterministic case the isochrones are 
singular objects which become smooth and well-defined 
due to noise; (iv) excitable systems which do not os- 
cillate without noise and therefore have no isophases, 
but the latter appear for the noise-induced dynamics. 



Furthermore, we have demonstrated applicability of the 
method to irregular experimental data. The definition of 
isophases in noisy systems has two potential application 
fields: in the data analysis, where it allows one to per- 
form a consistent phase reduction for signals with large 
amplitude variations, and in the synchronization theory, 
serving at a determining of phase responses to external 
kicks. 
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